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I]  abstract 
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A  seismic  location  study  has  been  conducted  in  a  continental  region,  Central  Alaska,  to  verify 
the  general  applicability  of  the  travel-time  anomaly  technique. 

First,  using  local  stations  and  estimating  event  depths  by  observing  pi),  a  set  of  earthquakes  4 
in  the  Central  Alaska  region  was  located  with  residual  travel  time  errors  attributable  to  reading 
error.  The  distribution  of  anomalies  for  teleseisraic  stations  indicates  that  the  Centra], Alaska 
region  is  composed  of  three  subregions  between  which  the  anomalies  vary  by  os  much  as  5 "seconds 
for  some  stations.  Some  events  on  the  edge  of  the  local  network  may,  however,  be  mislocatcd.  Such 
an  error  could  suggest  an  anomaly  change  where  there  was  none. 

The  size  of  one  of  the  subregions,  Alaska  Range,  is  about  ISO  km  by  ISO  km,  across  which 
constant  station  anomalies  are  apparently  valid  for  Ideating  using  a  combined  local  and  teleseismic 
net.  Within  this  subregion,  the  event  depths  vary  by  about  40  km,  from  85  km  to  125  km;  this  suggests 
that  no  significant  .anomaly  change  is  incurred  due  to  depth. 

Using  a  network  comprised  of  both  local  and  teleseismic  stations,  and  applying  as  corrections 
anomalies  determined  from  calibration  events  in  each  of  the  thrt.  subregions,  the  locations  of  the 
other  events  shift  less  than-S  km  from  the  local  network  epicenters  and  the  standard  deviations 
reduce  to  acceptably  low  values. 

A  comparison  of  locations  obtained  using  a  local  network  with  those  using  only  a  teleseismic 
network  car  be  reliably  obtained  for  only  two  events.  Tlie  shifts  from  the  local  locations  ar*  3  and 
11  km  with  anomalies,  as  contrasted  to  39  and  25  km  wit  -t  anomalies. 

There  is  an  indication  of  an  ave'rage  teleseismic  bias  of  13  km  to  the  northwast  over  most  of 
Central  Alaska,  thus  suggesting  a  large  anomaly  region.  Teleseismic  locations  obtained  with  a  set 
of  average  anomalies  agree  as  well  with  local  locations  as  do  locations  using  anomalies  from  events 
in  each  apparent  anomaly  sub-region.  3 
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ABSTRACT 

A  seismic  location  study  has  been  conducted  in  a  conti¬ 
nental  region,  Central  Alaska,  to  verify  the  general  applica¬ 
bility  of  the  travel-time  anomaly  technique. 

.  O'* 

First,  using  local  stations  and  estimating  event  depths 
by  observing  pP,  a  set  of  earthquakes  in  the  Central  Alaska 
region  was  located  with  residual  travel  time  errors  attribut¬ 
able  to  reading  error.  The  distribution  of  anomalies  for 
teleseismic  stations  indicates  that  the  Central  Alaska  region 
is  composed  of  three  subregions  between  which  the  anomalies 
vary  by  as  much  as  5  seconds  for  some  stations.  Some  events 
on  the  edge  of  the  local  network  may,  however,  be  mislccated. 
Such  an  error  could  suggest  an  anomaly  change  where  there  was 
none. 

The  si^e  of  one  of  the  subregions,  Alaska  Range,  is  about 
150  km  by  150  km,  across  which  constant  station  anomalies  are 
apparently  valid  for  locating  using  a  combined  local  and  tele- 
.seismic  net.  Within  this  subregion,  the  event  depths  vary  by 
about  40  km,  from  85  km  to  125  km;  this  suggests  that  no  signi¬ 
ficant  anomaly  change  is  incurred  due  to  depth. 

.  i  y  a 

Using  a  network  comprised  of  both  local  and  teleseismic 
stations,  and  applying  as  corrections  anomalies  determined  from 
calibration  events  in  each  of  the  three  subregions,  the  loca¬ 
tions  of  the  other  events  shift  less  than  5  km  from  the  local 
network  epicenters  and  the  standard  deviations  reduce  to 
acceptably  low  values. 

A  comparison  of  locations  obtained  using  a  local  network 
with  those  using  only  a  teleseismic  network  can  be  reliably 


obtained  for  only  two  events.  The  shifts  from  the  local  loca¬ 
tions  are  3  and  11  km  with  anomalies,  as  contrasted  to  39  and 
25  km  without  anomalies. 

There  is  an  indication  of  an  average  teleseismic  bias  of 
13  km  to  the  northwest  over  most  of  Central  Alaska,  thus 
suggesting  a  large  anomaly  region.  Teleseismic  locations 
obtained  with  a  set  of  average  anomalies  agree  as  well  with 
local  locations  as  do  locations  using  anomalies  from  events 
in  each  apparent  anomaly  sub-region. 
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INTRODUCTION 


Studies  of  travel-time  anomalies'  and  their  effects  on 
location  accuracy  and  network  stability  have  been  conducted 

•  •  B  . 

in  several  epicentral  areas.  For  example,  at  the  Neva-da  Test 

Site  , where  the  data  set  is  comprised  only  of  explosions,  the 

^ technique  of  using  travel-time  anomalies  as  corrections  can 

reduce  location  eWors  by  an  order  of  magnitude,  from  2  5  km 

to  2. 5-3.0  km;  furthermore,  location  Bias  can  be  shown  to  be 
1  »  -  f 
a  function  of"  the  particular. network'  used  for  locating 

(Chiburis,  1968;  Chiburis  and  Ahner,  1970).  }' 

p-  '  . 

I f  the  event  data  are  comprised  of  earthquakes,  location 

C  ,  *  9  f 

accuracy  cannot  be  assessed  in  an  absolute  sense  because  the 
true  locations  of  the  events  are  not  known..  That. a  significant 
■’  error  is  incurred  when  locating  without  t  ravel  ~"£iroe  anomalies  in 
a  region  suctf  as  the  Aleutian  Islands  has  been  demonstrated  by 
locating  LONG  SHOT,  an  underground  nuclear  explosion  on 
Amchitka  Island  in  the  Rat  Islands,  with  a  329  station  network. 
The  resultant  location  error  (with  depth  restrained)  exceeds3 
20  km  (Lambert  et  al.,  1970).  As  part  of  a  recent  study, 
Chiburis  and  Racine  (1971)  took  several  stable  subsets  of 
this  network  and  showe^i  that  thb  resultant  location  errors  of 
LONG  SHOT  were  variable  both  in  magnitude  and  direction, 
frequently  exceeding  30  km  (but  also  coincidentally  as  small 
as  1  km)  and  generally  random  in  directibh,  although  thjere 
was  an  apparent  northerly  bias.-  From  the/e  results,  it  was 
established  that  atcurate  locations  are  not  necessarily  ' 
achieved  simply  by  having  a  wjpll  distributed  network.  In  the  . 
same  report,  id  was  further  shown  that  location  consistency 
for  a  selected  set  of  earthquakes  across  the  AleutiarKlslands 
region  could  be  obtained  if  travel-time  anomalies  were 


« 


r  + 


■***&■* 


functionalized  for  each  station  recording  all  of  the  events; 
this  consistency  was  manifested  when,  by  using  network  subsets, 
a  small  cluster^of  locations  was  achieved  for  any  event,  and 
when  each  of  the  locations  displayed  acceptably  small  least- 
squares  .time  errors. 

¥'  t 

1  he  re fore ,  the  technique  of  applying  travel-time  anomalies 
even  if  the  anomalies  are  determined  from  £n  unknown/ or  mis- 


located  event,  is  valid  for  improving  location  accuracy 
(although  perhaps  only  relative  accuracy)  in  regions  such  as 
the  vNevada  Test  Site  and  the  Aleutian  Island?.  By/th  of  these 
regions  are  quite  different  geologically ;  thh~$/vada  Test 
Site  region  is(in  a  Basin  and  Range  province,  whereas  thb 
Aleutian  Islands  region  is  a  typical  island-arc  structure. 
Although  the  type  of  structure  in- the  vicinity  of  the  source  ( 
region,  should  have  no  effect  on  the  applicability  of  the 
anomaly  technique,  especially  if  only  teleseismic  stations  < 
are  used,  it  is  necessary  that  the  technique  be  tested  for 
events  occurring  in  a  ''continental -type"  s't-^uctu're  to  ^kiifther 
verify  the  procedures.  Therefore,  one  of  the  objectives  of  thp 
present  study  is  to  obtain  an  event  set  from  a  continental 
region,  such  as  Central  Xiaska,  and  investigate  .the  variability 
of  observed  travel-time  anomalies  at  both  local  and  teleseismic, 
stations  as  a  function  of  event  •position  within  th^e^region. 

Another  objective  is.  to  use.  events  whose  depths  are  known 
principally  frorf  identified  depth  phases  (pP)  and-  attempt  to 
discern  any  anomaly  instability  or  change  as  a  function  of 
depth.  #  '  - 

A  third  objective  is-  to  investigate  the  lo^tion  Capability* 
of  a  local  network  (all  stations  with  less  than  1*  000  km  epicen- 
tral  distances) ,  and  using  anomalies  determined  from  the  local 

network,  of  a  combined  local  and  teleseismic  network  and  of 
a  teleseismic  only  network. 


’  .  I 


REGION  DEFINITION  AND  DESCRIPTION  OF  THE  DATA 

defined  aTthat  ^  7'"*' '**“  Ce"tral  »«*«  is 

b  .  .  ,  3rea  between  Latitudes  62°N  and  66°N  and 

between  Longitudes  146°W  and  154°W  the  i 

Figure  1  show,  ♦  »  •  ™  shaded  portion  of 

g  l  shows  ti  region  of  interest,.  ^ 

The-  33  events  selected  for  analysis  were  of  magnitude  ' 

.  or  greater  as  reported  by  the  National  Ocean  Survey  (NOS) 
Preliminary  Determination  of  Epicenters  (PDE)  lists  for  the 
time  period  1964,  through  ISb^Figure  2  show  the  distri 
button  of  the  33  events  within'the  area  of  dnteres  al  \ 

Included  on  this  figure  are  the  event  name,\the  repo-ted  '' 
depths  and  7 ,000- foot  ..contours  ‘in  the  Alaska  Ran  e  ble  1 
gives  the  event  parameter  as  reported  by  NOS. 

-  thIhe  STi0^  SeleCted  for  ‘be  various  networks  total  i48f 
tTs?7;;e\SA  7  K°rld.Wide  Seismio  Statics  '• 

five  are  the  F  b"*  7'  Sels"lc ' Measurement  (LRSMKtations , 

the  Larv  a  ■  rl"8  a"d  A0  subarra>’  center  instrumentT'from 
VSLAWservatories  (Ofiev  7.  (USA)  ln  M°"tana-  fire 

rr.  rs-c  xr- 


*V> 


■**53K3Ktt  ja-s*  - -  *  - 

dna  eventJPS  shown  in  Table  ttt« 

station  +C  •  .  ’  he  reading  accuracy  by 

ype  is  estimated  to  be  as  follows-  LASA  and  Dh 
tory  ±0.1  sec;  LRSM  ±0.2  sec-  mss'1  +  n  c  ‘  7  Observa- 

sec,  mSS  ±0.5  sec;  U-  of  A  ±0.1  sec. 
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Figure  2.  Distribution  of  epicenters 


V  TABLE  I 


Event  Information  Reported  by  the  NationaL.Qcean  Survey  (NOS) 


Event 
Numbe  r 

Event 

Date 

Geographic 

Region 

Origin 

Time 

Latitude 

(NOS) 

Longitude 

(NOS) 

Magnitude 

(NOS) 

Depth 

(NOS) 

1 

29  Oct  68 

Tanana 

22 

16 ;  16 

65. 4N  ' 

150, 1W 

6.0 

7 

31  Oct  68 

Tanana 

00 

25:45 

65. 4N 

150. 1W 

4.5 

16 

3 

06  Feb  67 

Fai rbanks 

14 

48:40 

64. 8N 

147. 4W 

4.5 

53 

4 

21  Jun  67 

Fai rbanks 

18 

04:49.5 

64. 76N 

’  ,  147. 37W 

5.4 

17 

5 

21  Jun  67 

Fairbanks 

18 

13:02.9 

64 . 76N 

147. 37W 

5.6 

17 

6 

,23  Jun  67 

Fai rbanks 

11 

54:33,5 

64. 82N 

147. 45W 

4.6 

9 

7 

08  Sep  68 

Fai rbanks 

16 

22:58 

64. 8N 

147. 6W 

4.5 

12 

8 

07  May  63 

S.  of  Fairbanks 

05 

16:21 

64. 5N 

146.9W 

4.5 

33 

9 

30  Jun  63 

S.  of  Fairbanks 

05 

18:09 

64. 4N 

146. 9W 

4.5 

27 

10 

06  Apr  63 

Alaska  Range 

11 

19:23 

63. 4N 

149. 5W 

5.5 

39 

11  ■ 

06  Apr  63 

Alaska  Range 

12 

07:10 

63. 6N 

149. SW 

4.  5 

55 

12 

02  May  63 

Alaska  Range 

23 

13:13 

63. 2N 

148. 9W 

6.1 

81 

13 

22  Aug  63 

Alaska  Range 

03 

58:43 

63.  2N 

148. 5W 

4.6 

101 

14 

14  Dec  6  3 

Alaska  Range 

07 

51:08 

62. 7N 

149. 5W 

5.1 

95 

15 

1 4  May  6  4 

Alaska  Range 

11 

55:28 

62. 8N 

152. 3W 

4.6 

15 

16 

29  Jun  64 

Alaska  Range 

07 

21:33 

62. 7N 

152. 0W 

5.6 

33 

Q 

17 

06  Sep  64 

Alaska  Range 

"1 7 

36:44 

63. IN 

147. 7W 

4.8 

33 

18 

20  Nov  64 

Alaska  Range 

21 

27:'40 

63. 7N 

146. 5W 

4.6 

80 

19 

27  Nov  64 

Alaska  Range 

07 

47:08 

62. 6N 

151. 5W 

5.4 

113 

20 

21  Dec  64 

Alaska  Range 

18 

32:03 

63. IN 

150. 3W 

4.8 

111 

21 

08  Feb  65 

Alaska  Range 

03 

37:35 

63. 4N 

151. 7W 

4.5 

31 

22 

08  Mar  65 

Alaska  Range 

12 

04:21 

/ 

62. 5N 

150. 4W 

4.5 

104 

23 

26  Jun  65 

Alaska  Range 

23 

13^40 

62. 7N 

148. 7W 

4.8 

33 

24 

16  Oct  65 

Alaska  Range 

11 

4S:Wy 

63. IN 

150. 4W 

4.6 

84 

25 

24  Nov  65 

Alaska  Range 

08 

22:39 

63. 2N 

151. 0W 

5.1 

129 

26 

11  May  66 

Alaska  Range 

01 

26:24 

62. 8N 

150. 1W 

4.7 

91 

27 

31  Mar  67 

Alaska  Range 

04 
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The  principal  reasons  for  the  different  estimated  accuracies 
are  the  differences  in  seismogram  recording  speed;  however, 
overall  quality  of  the  seismic  traces  (thickness focus, 
etc.)  is  also  a  factor.  ' 


IMA 
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INITIAL  TRAVEL-TIME  ANOMALIES 


Using  the  NOS  reported  hypoc<?nter  parameters  listed  in 

Table  Ij  travel-time  anomalies  were  calculated  for, all  sta- 

v  tions  and  events.  Previous  experience  in  other  regions  has 

shown  that  the  anomalies  vary  significantly  with  location; 

therefore,  the  reported  locations  we r^  used  to*separate  the 

events  i.nto  three  .subregions  so  as  to  note  any  anomaly 

r  patterns'  across  Central  Alaska.  These  three  subregions  are 

-  referred  td  in  this  report  as  Tanana  (which  includes  2  events)  , 
.  P  r 

Fairbanks  (7  eVents),  and  Alaska  Range  ^,2  4  events)  ;  the 
dashe4" line  areas  in  Figure  2  show  the  space  relationship 
.  of  the  three  subregions. 

Travel-time  anomalies  are  computed  in  the  usual  manner; 
the  anomaly  at  station  i  relative  to  station  j  is 


Ai/i  -,Ti  -.tj  •  *  v4'-3? 


T-S 


where  T  is  station  arrival  time'  and  H  is  predicted  travel-time 
for  an  event  at  distance  A  and  depth  z  based  on  the  travel - time 
table  of  Herrin  (1968).  All  anomalies  in  this  report  are  made 
relative  to  st^tioh  UBSO  as  a 'matter  of  convenience.  The 
anomalies. by  event  for  several  stations  <are  given  in  Table  IV; 
these  stations  were  selected  to  illustrate  IfhcS  generally 
poor  consistency  observed.  With  anomaly  differences  within  a 
subregion  larger  than  4  seconds  at  teleseismic  stations  and 
6  seconds  at  local  stations,  seismogram  reading  errors  can 
immediately  be  eliminated  as  the  sole  source  of  the  discrepancies. 
Previous  experience  has  shown  that  the  chief  causes  of  anomaly 
Inconsistency  are  epicenter  mislocations  and  depth  errors 
(Chiburis  and  Dean,  1965  ;  Chiburis,  1968^Chiburis  and  Ahner,  1970; 

...  .  v 
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TABLE  IV 


•  ^ 


Travel -Time 

Anomalies 

at  Selected 

Stations 

As 

Computed  from 

NOS  Hypocenters  Listed  in  Table  I 

Event 

No. 

Date 

COL 

NP-NT 

SCP 

TNN 

1 

29  Oct  68 

-2;  44 

-5,60 

-3.04 

-0.17 

2 

31  Oct  68 

-2.83 

-6.18 

-2.41 

-1.51 

Ave  rage 

-2.64 

-5.89 

-2.73 

-0.84 

Standard  Deviation 

0.28 

0.41 

0.45 

0.95 

3 

06  Feb  67 

-3.18 

-1.34 

-- 

+  1.62 

4 

21  Jun  67 

-5.55 

-2.14 

- 

-4.  56 

5 

6 

21  Jun  67 

2  3  Jun  6  7 

-2.70 

-3.09 

-3.23 

7 

08  Sep  68 

-- 

-2.76 

-- 

-2.20 

8 

07  May  63 

-- 

-- 

-- 

-- 

9 

30  Jun  63 

• 

-- 

-0.63 

-- 

Average 

-3.81 

-2.33 

-0.63' 

-2.0^ 

Standard  Deviation 

1.53 

0.77 

-- 

2.66 

10 

06  Apr  63 

-- 

-- 

-1.28 

11 

06  Apr  63 

-- 

-- 

-1.05 

12 

02  May  63 

-- 

-- 

-0.26 

-- 

13 

22  Aug  63 

-- 

+  1.90 

-0.55 

-- 

14 

14  Dec  63 

-- 

-1.50 

-3.44 

-- 

15 

14  May  64 

-0.17 

+2.17 

-0.34 

-- 

16 

29  Jun  64 

-- 

-1.44 

-1.25 

-- 

17 

06  Sep  64 

-- 

-1.29 

-- 

18  ' 

20  Noy  64 

-2.27 

+0.51 

-0.25 

-- 

19 

27  Nq 4  64 

-0.76 

-1.17 

-0.58 

-- 

20 

21  Dec  64 

+0.14 

+  1.31 

-- 

-- 

21 

08  Feb  65 

-2.75 

-0.99 

-4.64 

-- 

22 

08  Mar  65 

-1.68 

-- 

-  - 

23 

26  Jun  65 

+  1.96 

+  2.21 

-1.42 

-- 

24 

16  Oct  65 

-1.75 

-0.17 

-1.69 

-- 

25 

24  Nov  65 

-0.91 

+0.68 

-1.37 

-- 

26 

11  May  66 

-1.41 

+0.71 

-1.25 

-- 

27 

31  Mar  67 

-- 

-1.51 

-- 

-1.05 

28 

05  May  67 

-1.78 

-1.71 

-1.45 

-0.70 

29 

06  Jul  67 

+  0.41 

-- 

-- 

-- 

30  11  Oct  67 

-1.11 

-  -- 

-1.05 

-0.46 

31 

32’  ^ 

10  Nov  b7 

04  De  c  ^7 

-1.  35 

-1.13  > 

-1.13 

-0.77 

-0.22 

-0.22 

-o.-rf'" 

+0.24 

33 

28  Dec  68 

-3.41, 

-2.60 

-1.90 

Average 

-1.  07 

-0.24 

-1.32 

-0.77 

Standard  Deviation 

1.31 

1.47 

1.10 

0.70 

Cve  rail 

Ave  rage 

-1.61 

-1.06 

-1.42 

-1.22 

Standard  Deviation 

1.61 

2.13 

1.12 

1.63 

) 


Chiburis  and  Racine,  1971).  To  eliminate  possible  errors  in 
the  reported  depths,  each  of  the  events  was  analyzed  by 
aligning  the  P  waves  at  all  stations  as  a  function  of  distance 
and  then  selecting  apparent  pP  phases  based  on  coherent  energy 
across  the  network.  The  assistance  which  such  a  technique  can 

provide  is  well  known  and  was  demonstrated  by  Chiburis  and 

.  s, 

Ahner  (1969).  The  observed  (pP-P)  time  differences  were 
converted  to  depth  using  the  relation  of  Jeffreys  and 
Shimshoni  (1964).  A  comparison  between  the  depths  thus 
obtained  and  the  NOS  reported  depths  is  shown  in  Table  V. 

Of  the  33  event  depths,  twenty  disagree  by  less  than  7  km 
and  only  seven  disagree  by  more  than  20  km  of  whi^ch  four 
disagree  by  more  than  60  km. 
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TABLE  V 

fl 


A 

Comparison  Between  TIOS 

Reported  and 

Depths 

Determined  by  pP 

Observations 

Eve  n't 

NOS  Depth 

pP  Depth 

<Jz 

No. 

Date 

(km) 

(km) 

(km) 

1 

29 

Oct 

68 

7 

13 

+6 

2 

31 

Oct 

68 

16 

16 

0 

3 

06 

Feb 

67 

S3 

28 

-25 

4 

21 

Jun 

67 

17 

16 

-1 

S 

21 

Jun 

67 

17 

13 

-4 

6 

23 

Jun 

67 

9 

16 

+  7 

7 

08 

Sep 

68 

12 

8 

-4 

8 

07 

May 

63 

33 

13 

-20 

9 

30 

Jun 

63 

27 

9 

-18 

10 

06 

Apr 

63 

39 

108 

+  69 

11 

06 

Apr 

63 

-55 

115 

+  60 

12 

02 

May 

63 

81 

85 

+  4 

13 

22 

Aug 

63 

•  101  . 

78 

-23 

14 

14 

Dec 

63 

95 

108 

+  13 

15 

14 

May 

64 

IS 

16 

+  1 

16 

29 

Jun 

64 

33 

"  21 

-12 

17 

.06 

Sep 

64 

33 

134 

+  101 

18 

20 

Nov 

64 

80 

85 

+  5 

19 

27 

Nov 

64 

113 

113 

0 

20 

21 

Dec 

64 

111 

108 

-3 

21 

08 

Feb 

65 

31 

15 

-16 

22 

08 

Mar 

65 

104 

105 

+  1 

23 

26 

Jun 

65 

33  • 

134 

+  101 

24 

16 

Oct 

65 

84 

97 

+  13 

25 

24 

Nov 

65 

129 

134 

+  5 

26 

11 

Nov 

66 

91 

85 

-6 

27' 

31 

Mar 

67 

83 

82 

-1 

28 

05 

May 

67 

102 

100 

-2 

29 

06 

Jul 

67 

59 

54 

-5 

30 

11 

Oct 

67 

115 

113 

-2 

31 

10 

Nov 

67 

90 

93 

+  3 

32 

04 

Dec 

67 

96 

108 

+  12 

33 

28 

De  c 

68 

80 

82 

+  2 

(  ' 
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RELOCATIONS  USING  A  LOCAL  NETWORK 

If  the  pP-determined  depths  are  accepted  as  being  correct, 
a  local  network  can  be  used  to  relocate  the  epicenters  of  the 
12  events  satisfactorily  recorded  by  it,  restraining  the  depths 
to  the  accepted  values  and  applying  no  anomalies  as  corrections. 

The  geometry  of  the  local  network  used  for  this  purpose  is 
shown  in  Figure  3. 

The  location-shifts  from  the  12  NOS  reported  epicenters  and 
the  standard  deviations  of  the  least-squares  time  fits  for  the 
local  network  are  given  in  Table  VI  ;  also  indicated  are  the 
stations  used  in  making  each  of  the  locations.  We  used  the 
Herrin  (1968)  travel  time  curves.  A  significant  result  in  this 
table  is  the  improvement  achieved  in  the  standard  deviation  of 
the  time-fit  location  when  compared  with  the  standard  deviation 
determined  from  the  input  location.  Except  for  the  Tanana 
subregion,  where  the  improvement  is  negligible  and  where  the 
local  network  control  is  poorest,  the  resultant  standard 
deviations  approach  those  which  one  would  expect  from  reading 
errors  alone. 

Some  -ither  points  relevant  to  Table  VI  :  the  NOS  reported  loca¬ 
tion  of  the  06  February  1967  Fairbanks  event  was  clearly  in  error. 
The  standard  deviation  for  six  local  stations  is  nearly  4  seconds 
as  computed  from  the  NOS  epicenter  and  pP  depth.  Since  the  NOS 
location  was  obtained  with  a  depth  of  53  km  instead  of  28  km, 
it  is  probable  that  the  depth  error  was  responsible  for  most 
of  the  mislocation.  Using  the  local  network,  the  solution 
obtained  for  this  event  with  a  restrained  depth  of  28  km  is 
now  consistant  with  the  other  events  in  the  Fairbanks  subregion. 
Excluding  the  06  February  1967  event,  the  average  shift  of  the  event 
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TABLE  VI 

Location  Results  Using  a  Local  Network 
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set  is  about  13  km  from  the  NOS  epicenters.  A'  geographic  pi  ,t 
o  the  location  shift  vectors  is  shown  in  Figure  4.  except  or 
two  events  in  the  Fairbanks  subregion,  the  directions  of  th . 
location  shifts  from  the  NOS  locations  vary  between  east  and 
southeast.  This  apparent  bias  cannot  be  explained  in  terms  of 
network  instability  or  distribtuion.  The  relation  between 
network  geometry  and  location  shift  is  shown  in  Figure  5  in 
w  ich  the  station  distances  are  plotted  as  logarithmic  vectors 
along  their  azimuth  lines  and  the  location  shifts  as  lihear 
vectors.  The  results  in  this  figure  do  not  indicate  any 
correlation  between  the  direction  of  event  shift  and  the 
distribution  of  the  network  used  for  locating  that  event.  In 
this  regard,  however,  it  was  shown  in  an  earlier  report  '(Chibur 
and  Ahner,  1970)  that  location  bias  is  a  function  of  the 
particular  network  defined  for  locating,,  and  that,  if  a 
-™Stant  net“°rk  is  used  to  locate  a  set  of  events,  the 
direction  and  magnitude  of  the  shifts  for  all  events  within 
an  anomaly  region  are  closely  similar.  It  was  further  shown 
that,  by  carefully  selecting'*  network,  virtually  any  bias 
direction  and  magnitude  one  desires  could  be  obtained  for  the 
same  event  set.  Although  a  constant  network  is  not  used  in  the 
present  case,  this  "network  bias"  could  still  be  a  significant 
-factor  for  producing  the  observed  southeasterly  shifts 
because  the  12  events  were  located  using  networks  composed 
largely  of  the  same  stations.  In  particular,  stations  Pjti,  TNN 
CM  an  WH2YK  were  used  in  all  12  of  the  locations;  statins 
COL,  BIG  and  BLR  were  used  in  nine  of  the  locations.  If  the 
travel-time  anomalies  for  these  stations  happen  to  be  large 
the  similarity  among  the  'spveral  networks  using  them  may  be’'’ 
enough  to  explain  the  observed  bias.  .-  • 

In  another  report  (Ghiburis  and  Racine,  1971),  it  was 
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shown  that,  although  network  bias  usually  predominates,  a 
\  secondary  and  sometimes  consistent' bias  can  exist  almost' 
regardless  of  the  network  used,  at  least  f  qr-'Jexpl_os ions 
in  the ’Aleutian  Islands,  and  the  southeasterly^fiif ts 
observed  on  Figures  4  and  5  may  possibly  be  due  to  this 
effect.  In  particular  the  local  network-  locations  may  be 
correct  and  the  NOS  epicenters  biased  to  the  northwest. 

This' would  be  'consistent  with  stf\?  northward  shift  of 

LONG-SHOT  observed  with  a  telese&smic  network , ‘and  is“the 

-  *  5 

Explanation  we  prefer.  An  implication  would  be  that  the  > 

'central^  Alaskan  bias  can  be  traced  fo  the  Aleutian  under-  .  ^ 

-  f  4  ’  „ 

thrusting  lithospheric  $lab.  However,  the  Central  Alaskan 
region  is  a  structurally  complex  area  and  the  bias  may  In 
fact  be  related  to  spme  other  ^structural  ‘feature. 

It  is  important  at  this  point  to  compare  the  residQals 
and  least-squares  time  fits  obtained  from  the  local  network 
solutions  to  those  obtained  from  the  NOS  solutions.  Using- 

v 

the  N5S  hypocenters,  a  residual  at  station  i  for  "every;  j  can 
be  calculated  ^is 


rJ  =  tI  -  T^  -  H-?  (A ,  z) 

i  n  c\  i  '  ' 


where  T  is  the  event  origin  time  and  the  other  terms  ate  as 
o  ~  <  ■ 

defined  before.  The  robt-mean-square'  of  the  residuals  for  the 
total'  event  and  local  station  set  is: 

L  j-1  i=l  L  j=l*  ■  J  . 
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where  N.  is  the  number  of  local  stations  recording  the  ith, 

J  -  ,  \  ^ 

event.  The  value  of  a NOS  ^*^58  sec  for 'these  data.  Assuming 

that  this  figure  generally  reflects-  the  goo<dness-of- f  it  to 

the  travel-time  model  used  (in  this^ c^se /.Herrin  (1968)),  it 

can  be  stated  th\t  the  fit.  is  poor.  However,  since  many  of  f 


the  stations  within  the  local  network  shown  in  Figurtev3  were 
not  included  in  the  NOS  solution  for  these  events,  an  "origin¬ 
time  bias"  for  each  event  should  be  remcWell  from  each  re3i  ’ual 
before  combining  the  set;  that  is 


=  T-l  -  TJ  -  H?  -  RJ.  ’ 

1X01  * 


where  ^R-*  is  the  mean  |p f  the  local-ststion  residuals  for  the  j  th 
event:-  / 

'12  *  RJ 

— r-  5—*  1  ( 

rj  -  2_,  ~  * '  y  » ' 

i  - 1  N  ' 


The  standard  deviation  is  then  calculated  as 
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The  value  of'< 


^NOi/  :*'s  ^  sec,  a  reduction  of  0.48  sec  compared 


to  the  straight  root-mean-square .  This  is  still  high  enough  to 
indicate  a  poor  fit. 


If  now  the  standard  deviation  of  residuals  is  calculated 
on  the  basis  of  solutions  obtained  using  only  the  local  network, 


the  value  of  a 


LOC 


0.511  sec.  This  result  suggests  that  the 
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crustal  and  upper  mantle  model  used  for  the  Herrin  (1968) 
tables  is  a  reasonably  good  average  representation  of  the 
crustal  and  upper  mantle  structure  in  Central  Alaska;  this 
suggestion  is  in  essential  agreement  with  conclusions  drawn 
by  Berg  et  al..(1967),  who  concluded  that  'the  P-wave  velocity 
in  the  upper  part  of  the  crust  in  the  Fairbanks  area  is  about 
5.56  km/sec,  which  is  also  near  the  value  of  6.0  km/sec  used 
in  -the  Herrin  (1968)  model. 

To  test  further  the  appropriateness  of  the  crustal  model, 
the  local  network  can  be  used  to  locate  the  events  with  the 
depth  parameter  included  in  the  least-squares  adjustment.  The 
results  of  doing  this  are  given  In  Table  VII.  The  Tnaj  ority  of  ^ 
events  move  only  a  few  kilometers  in  depth  compared  -to  the 
locations  obtained  with  the  depths  restrained  to  a  value  based 
on  pP.  The  exceptions  to  this  depth  stability  are  the  two 
Tanana  events,  which  change  in  depth  by  "2  km  and  25  km  upward 
and  in  epicenter  both  by  about  30  km  northward;  the  reason 
for  this  is  that  these /events  are  not  within  the  geometric 
confines  of  the  local  network  (Figure  3),  so  the  soJutions 
tend  to  move  from  the  network  in  Ibe  direction  of  least 
control.  Two  otheb  possible  exceptions  to  the  depth  stability 
are  the  events  31  March  1967  and  28  December  1968  in  the  v 
Alaska  Range.  The  unrestrained  solution  for  31  'March  1967  , 

changes-' by  31'kmyupward  in  depth  and  12  km  east  in  epicenter; 
however,  this  could  be  e'xplainei  by -the  low  number  of  recording 
stations  whi<:h  provides  only  a  single  degree-  J>f  freedom  in 
the  least-squards  procedure.  The  28  De cemberv^968  event  shifts  , 

18  km  upward  in  depth  and  13  km  east-southeast  in  epicenter; 

’  i-t  is  riot  entirely  clear  why  this  event  shifts,  as  the  network 
is  well  distributed  about  the  epicenter Qdnd^ the  standard  \  £ 

,  „  *  r. 
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deviation  for  the  restrained  case  (ap  is  small.  Conceivably, 
this  event  may  actually  be  shallower  than  the  depth  indicated 
by  pP;  if  pP  was  mispickpd  by  about  3  sec,  the  depth  would 
be  in  the  general  range  of  the-  depth-free  solution.  This  possi¬ 
bility  will  be  further  discussed  in  the  next  section. 

1  4* 

If  the  12  restrained  hypoceitters  determined  by  the  local 

/ 

network  are  now  used  to  compute  the  anomalies,  the  results 

show  better  consistency  at  most  stations.  For  example,,  a 

comparison  of  the  anomalies  computed  from  the  NOS  location  to 

those  obtained  from  the  local  network  is  g'iven  in  Table  VIII 

for  stations  TNN,  COL  and  WH2YK,  which  are  part  of  the  local 

network,  and  for  stations  NP-NT  and  OXF,  which  are  teleseismic 

to  Central  Alaska  and  recorded  a  substantial  number  of  the 

events.  The  improvement  in  the  consistency  of  the  local  stations 

for  the  Fairbanks  subregion  is  expected  because  of  the  mislocated 

event  06  February  1967.  Because  small  depth  errors  are  not 

significant  for  teleseismic  stations,  the  improvement  in  anomaly 

consistency  at  NP-NT  and^OXF  is  less  significant  than  at  the 

local  stations.  The  overall  consistency  suggests  that  the  event 

locations  within  a  subregion  determined  with  a  local  networks 

/  ^  , 

are  more  accurate  (at  least  relative  to  one  another)  than  the 
reported  locations.  • 

One  could  not  reject  the  hypothesis  that  all  of  Central 
Alaska  is  one  anomaly  region  from  the  OXF  data.  However  for 
NP-NT.  (the  closest  teleseismic  station)  the  Tanana  anomalies 
(from  the  events  with  the  worst  local  network  control)  clearly 
define  a  different  region.  Overall,  these  data  leave  open  the 
possibility  that  all  of  Central  Alaska  is  a  single  anomaly  region 
with  respect  to  teleseismic  stations  farther  away  than  NP-NT. 

This  possibility  has,  of  course,  already  been  suggested  by  the 

•p 

consistent  northwest  bias  of  the  NOS  location  with  respect  to 
the  local  locations. 


TABLE  VIII 

Comparison  of  NOS  and  local-network  anomalies  for  Selected  Stations 
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TELESEISMIC  LOCATIONS  USING  CALIBRATION  EVENTS 

* 

to. 

If  a  calibration  event  is  selected  from  each  of  the  three 
subregions  in  Central  Alaska  and  anomalies  at  both  local  and 
teleseismic  stations  are  computed  from  the  restrained  hypo- 
centers  determined  in  the  previous  section,  and  if  the 
anomalies  are  used  as  corrections,  the  events  in  the  respective 
subregions  can  be  located  relative  to  the  calibration  events. 

In  this  way,  anomaly  stability  can  readily  be  evaluated  as  a 
function  of  distance  between  subregions  and,  to  a  lesser 
*  extent,  as  a  function  of  dep  th. 

The  events  selected  for  calibration  purposes  are  21  June 
1967  Fairbanks  (to  locate  06  February  1967,  23  June  1967  and 
08  September  1968);  29  October  1968  Tanana  (to  locate  31  Octobe 
1968;  and  10  November  1967  Alaska  Range  (to  locate  04  December 
1967  j  11  October  1967,  31  March  1967  ,  28  December  1968  and 
05  May  1967).  The  location  results  of  subregion  calibration 
using  confined  local  and  teleseismic  stations  are  given  in 
Table  IX  for  both  restrained  and  unrestrained  solutions.  In 
the  Fairbanks  subregion,  the  three  events  shift  about  5  km 
in  random  directions  and  the  standard  deviations  of  solution 
(a  )  reduce  to  about  0.7  sec,  less  than  half  the  standard 
deviations  obtained  without  anomalies  (oQ) .  The  depths  do  not 
Change  significantly  when  the  solutions  are  run  depth  free. 

In  the  Tanana  subregion,  the  single  event  shifts  about  7  km 
and  the  depth  changes  by  2  km.  Significantly,  the  standard 
deviation  is  now  reduced  to  about  0.5  sec.  When  compared  to 
the  standard  deviation  of  1.380  sec  from  the  restrained  hyj>o- 
center-  determined  by  the  local  network  (oq)  ,  it  is  clear  that 
the  anomalies  from  the  calibration  event  in  this  area  are 


TABLE  IX 

Subregion  Calibration  Results  for  Central  Alaska 
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correcting,  at  least  to  the  first  order,  the  imperfect  earth 
model  relevant  to  the  Tanana  subregion. 

For  the  Alaska  Range  subregion,  when  using  the  10  November 
1967  event  for  calibration  purposes,  the  average  restrained 
shift  of  the  five  events  is  only  about  2.5  km  in  a  random 
direction  and  the  standard  deviations  reduce  to  below  0.6  sec. 
fthen  the  solutions  are  allowed  to  run  depth  free,  the  average 
shift  is  less  than  2  km  and  the  standard  deviations  are  all 
below  0.4  sec  except  for  05  May  1967  which  is  farthest  from 
the  calibration  event.  Since  the  size  of  the  Alaska,  Range 
subregion  is  about  150*  km  by  150  km,  certainly, larger  than  the 
70  km  by  25  km  accepted  in  earlier  studies  (Chiouris  and  Ahner, 
19703  ,  it  is  expected  that  the  05  May  1967  event  would  have  the 
poorest  fit.  Foi  the  unrestrained  solutions  in  the  Alaska  Range, 
the  depths  of  the  events  change  only  a  small  amount  (less  than 
5  km),  except  for  two  events,  31  March  1967  and  11  October  1967, 
which  both  go  deeper  by  9|  km.  The  standard  deviations  of  these 
events  are  reduced  from  0.585  to  0.398  sec  and  0.517  to  0.274  sec 
respectively.  It  may  be  that  the  depth-free  solutions  are  yield¬ 
ing  more  reliable  depths  for  these  two  events,  although  it  may 
be  that  the  anomalies  are  changing  slightly  across  the  region. 

It  is  impossible  to  tell  if  these  two  events  are  correctly 
located  in  three  dimensions,  because  the  standard  deviations 
may  reduce  only  fortuitously  when  using  10  November  lv967  as  the 
calibration  event.  For  example,  if  the  31  March  1967^vent, 
restrained  to  a  depth  of  82  km,  is  used  to  calibrate  the 
11  October  1967  and  28  December  1968  events,  the  epicenters 
shift  respectively  about  1  km  and  14  km  with  depth  shifts, of 
4  km  downward  and  6  km  upward  respectively;  the  standard 
deviations  reduce  to  0.268  and  0.019  sec,  an  almost  perfect 
time  fit  for  the  latter.  Only  further  studies  will  remove 
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the  ambiguity. 


If  the  local  stations  are ' removed  from  the  networks  and 

only  teleseismic  stations  are  used,  the  results  of  locating 
/ 

both  with  and  without  anomalies  are  as  given  in  Table  X.  In 

9 

the  Fairbanks  subregion  the  shifts  are  large,  but  this  can  be 

9  , 

explained  by  poor  teleseismic  network  distribution.  Although 
the  networks  have  a  large  aperture  if  only  extreme  azimuths 
are  considered  (Network  Aperture  in  Table  X),  the  aperture 
is  much  less  if  only  one  or  two  stations  are  excluded.  A 
network  of  this  type  is  known  to  give  poor  locations.  For 
example,  the  01  February  1967  event  was  located  with  a  network 
having  an  aperture  of  107°;  if  one  station  is  deleted,  the 
aperture  (n-1)  reduces  to  42°  and  if  two  stations  -are  deleted, 
the  aperture  (n-2)  is  24°.  .  • 

In  the  Tanana  subregion,  the  teleseismic  network  used  to 
locate  the  31  October  1968  event  is  fairly  stable  down  to  the 
(n-2)  aperture  and  the  location  shift  is  about  10  km.  However, 
the  local  network  for  thi-s*  event  is  poor  with  apertures  of  157°, 
73°  and  55°  down  to  n-2  stations.  In  the  Alaska  Range  subregion, 
only  two  of  the  events,  04  December  1967  and  05  May  1967  have 
well-distributed  teleseismic  networks.  The  shifts  are  respectively 
3.2  km  and  11.1  km;  however,  as  pointed  out  earlier,  the  05  May 
1967  event  is  farthest  from  the  calibration  event  and  is  not 
expected  to  be  as  well -located.  Although  the  events  11  October 

M 

1967  and  31  March  1967  have  shifts  of  about  6  km  and  9  km 
respectively,  this  may  be  coincidental  because  the  (n-2)  tele¬ 
seismic  apertures  are  41°  and  20°,  hardly  large  enough  to  be 
classed  as  well-distributed.  The  28  December  1968  event  has  only 

’  f  ,  T 1 

three  stations  for  locating  (no  degrees  of  freedom  for  a 
restrained  location)  and  displays  the  largest  shift. 
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or  OS  May  1967  which  is  0.885  sec  If  w  )  .  • 

well-distributed  local  net  for  th  '•  ■  °  ret|Uire  a  . 

discard  the  10  km  shift  o'  the  T  lnltlal  locatlon  -we  will 
average  ertof  of  7  km.  e''e"t  multiai  ’»?.  a" 

When  we  compare  the  shifts  of  these  two  Alaska  R* 
obtained  without  anomalies,  the  shifts  in  i  E°  eVent5 

*.  3  km  and  from  15  km  L  Z  *?*■“  ^ 

the  shift  reduces  from  56  km  to  10  km!  ™  eVent* 

Also  in  Figure  5  we  have  plotted  the  sh.'ftc 

central  Alaskan  events.  We  see  that  rm  t-v, 
average  the  displacement  from  the  local  ,V,  ?  '• 

to  that  obtained  by  a  different  ca  b  ^  “  e«“al 

of  9.7  km.  ■  .  respecjtvfly.  fop  an  average  .  ' 

.  -  ‘  «* 

Thus,  although\he  reliable  results  are  sparse-  atpl 
location  error  of  ^-n  u  parse,,  a  teleseismic 

0t  3  11  km»  when  anomalies  are.  applied  i  e  *, 
expected  result  a«  ,  ^pried,  is  the 

h  a  result,  as  compared  to  an  iverage  error  nf  «  v  .  ,  * 

anomalies  for  two  Alaska  d„  8  >  *  of  32  km  without 

40  km  if  t  events  and  an  average  error  of 

40  km  if  the  Tanana  event  is  included. 
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CONCLUSIONS 


From  a  study  of  earthquakes  occurring  in  Central  Alaska, 
a  Continental  region  as  .opposed  to  an  island-arc  region  pr  a 
•"bas  in- and- range  region,  it  has  beeb  demonstrated  that  the  anomaly 
technique  is  valid,  especially  if  the  region  is  divided *into 
subregions  according  to  their  anomaly  stability. 

It  was  show)  that  a  local  network  could  ach ie ve • cons i s t^n t 
location  solution. »  as  regards  standard  deviation,  or  goodness- 
of-fit,  relative  to  the  Herrin  (1968)  travel-time  tables.  From 
these  local  locations,  anomalies  could  be  determined  for  both 
local  and  teleseismic  stations  and  applied  to  obtain  consistent 
so luiions  of  other  pvents  in  the  subregions. 

From  the  geographic  distribution  ot  the  events,  there  are 
three  anomaly  subregions  in  Central  Maska  (Tanana,  Fairbanks'  v 
and  Alaska  Range)  between  which  the  anomalies  change  significantly 
but  within  which  the  anomalies,  when  applied  as  corrections,  are*** 
apparently  correcting  the  inadequate  earth  model  (travel-  ime 
table).  There  i^s.  some  possibility  that  some  of  the  separation 
into  teleseismic  subregions  is  forced  by  poor  locations  in  the 
Tanana  region.  I  he  size  of  the  Alaska  Range  subregion  is  large,  . 
approximately  150  km^by  150  km,  where  five  events,  calibrated 
by  a  sixth,  yield  consistent  locations,  both  in  epieente-^  ajrd 
depth. 

There  are  some  suggestions  i'e.^U)  a  consistent  13  km 
southeast  bias  with  respect  to  NOS  epicenters  and  (2)  fair 

teleseismic  locations  resulting  from  averaged  anomalies,  tliat  -r 

# 

the  teleseismic  anomaly  region  even  larger  than  the  Alaska 
Range  subregion.  "  '  v 


The  location  errors,  obtained  by  applying  sub- region 
anomalies ,  are  less  than  5  km  based  on  nine  events  when  using 
local /and  teleseisrfi c  stations,  and  are  3-11  km  basetf  on  three 
event/s  when  losing  only  teleseismic  stations.  \ 
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